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We propose a way for generating n-qubit Greenberger-Horne-Zeilinger (GHZ) entangled states 
with a three-level qubit system and (n — 1) four-level qubit systems in a cavity. This proposal 
does not require identical qubit-cavity coupling constants, and thus is tolerant to qubit-system 
, parameter nonuniformity and nonexact placement of qubits in a cavity. The proposal does not 

. require adjustment of the qubit-system level spacings during the entire operation. Moreover, it is 

shown that entanglement can be deterministically generated using this method and the operation 
time is independent of the number of qubits. The present proposal is quite general, which can be 
applied to physical systems such as various types of superconducting devices coupled to a resonator 
or atoms trapped in a cavity. 
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Quantum entanglement plays a crucial role in quantum information processing. Entanglement of ten photons [1], 
eight ions [2], three spins [3], two atoms in microwave cavity QED [4], or two excitons in a single quantum dot 
[5] has been demonstrated in experiments. In addition, experimental preparation of three-qubit entanglement with 
^h. superconducting qubits or a superconducting qubit coupled to two microscopic two-level systems has been reported 
recently [6-8]. Although multi-particle entanglement was experimentally created in photons and trapped ions, it is 
still greatly challenging to create multi-qubit entanglement in other important physical systems. 

As is well known, multi-qubit GHZ (Greenberger-Horne-Zeilinger) entangled states are of great interest in the 
foundations of quantum mechanics and measurement theory, and significant in quantum information processing [9], 
■ quantum communication [10-12], error correction protocols [13], and high-precision spectroscopy [14]. Over the past 
ten years, based on cavity QED technique, many different theoretical methods for creating multi-qubit GHZ entangled 
states with atoms and superconducting qubits have been presented [15-25]. For instances, (i) the proposals in [15-17] 
for implementing a GHZ entangled state are based on identical qubit-cavity coupling constants; (ii) the approaches 
in [18-20] are aimed at probabilistic generation of a GHZ state; (hi) the method presented in [17] is based on the 
I ■ use of an auxiliary qubit, measurement on the qubit states, and adjustment of the qubit level spacings during the 
t-H I entire operation; (iv) the proposals in [21,22] are based on photon detection outside the cavity; and (v) the proposals 
in [23-25] are based on sending qubits (i.e., atoms) through a cavity. These proposals are important because they 
opened new avenues for creating multi-particle entanglement. 

In this paper, we will focus on a situation that the qubit-cavity coupling constants are nonidentical. This situation 
often exists in superconducting qubits coupled to a cavity or a resonator. For solid-state devices, the device parameter 
nonuniformity is often a problem, which results in nonidentical qubit-cavity coupling constants though qubits are 
placed at locations of a cavity where the magnetic fields or the electric fields of the cavity mode are the same. In 
addition, the nonidentical qubit-cavity coupling constants may also result from nonexact placement of qubits in a 
cavity. In the following, our goal is wish to present a way for deterministic preparation of a n-qubit GHZ state 
with tolerance to the qubit-system parameter nonuniformity and nonexact placement of qubits in a cavity. As shown 
below, this proposal also has these advantages: (i) The operation time is independent of the number of qubits in 
the cavity; (ii) When compared with the approach in [17], no auxiliary qubits, no measurement on the qubit states, 
and no adjustment of the qubit level spacings during the entire operation is needed (note that adjustment of the 
level spacings of the qubits during the operation is not desired in experiments and may cause extra errors); (iii) No 
adjustment of the cavity mode frequency is required during the entire operation; and (iv) There is no need of photon 
detection. This proposal is quite general, which can be applied to various types of superconducting qubits and atoms 
trapped in a cavity. 

This paper is organized as follows. In Sec. II, we briefly review the basic theory of a three-level quantum system 
or four-level quantum systems coupled to a single-mode cavity and/or driven by classical pulses. In Sec. Ill, we show 
how to generate an n-qubit GHZ state with one three-level quantum system and (n — 1) four- level quantum systems 
in a cavity. In Sec. IV, we compare our proposal with the previous ones. In Sec. V, we give a brief discussion of the 



2 



Q, 



(a) 



(b) 



(c) 



5: 



|2) 









#7 



A C7 



(d) 

FIG. 1: (Color online) (a) and (c) System-pulse resonant interaction for qubit system 1. In (a), the pulse is resonant with the 
|1) <-¥ 2) transition; while in (c) the pulse is resonant with the |0) o |1) transition, (b) System-cavity resonant interaction 
for qubit system 1. The pulse is resonant with the |1) |2) transition, (d) System-cavity-pulse off-resonant Raman coupling 
for n — 1 qubit systems (2,3, n) in a cavity. For simplicity, we here only draw a figure for qubit system j interacting with 
the cavity mode and a classical pulse (j = 2,3, ...,n). 5j = A c j — Aj is the detuning of the cavity mode with the pulse, 
Aj = ui 3 31 — ujj is the detuning between the pulse frequency uij and the |1) -R- |3) transition frequency oj^i of qubit system j. 
Note that gj, 5j, A c j, and Aj may be different for qubit systems (2,3, ...,n) due to nonidentical level spacings of the qubit 
systems which are caused by the nonuniformity of the qubit system parameters. The coupling constant gj may also vary with 
qubits due to nonexact placement of qubits in a cavity. The Rabi frequency of the pulse applied to qubit system j is denoted 
by Qj . 



experimental issues and possible experimental implementation with superconducting qubits coupled to a resonator. 
A concluding summary is given in Sec. VI. 

II. BASIC THEORY 

In this section, we will introduce three types of interaction of qubit systems with the cavity mode and/or the pulse. 
The results presented below will be employed for generation of a multi-qubit GHZ state discussed in next section. 

A. System-pulse resonant interaction 

Consider a three-level qubit system (say , qubit system 1) driven by a classical pulse. Suppose that the pulse is 
resonant with the transition |1) O |2) of the qubit system 1 but decoupled from the transition between any two other 
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levels [Fig. 1(a)]. The interaction Hamiltonian in the interaction picture is given by 

frj = ft(fi r e i *|l)<2|+H.c.) ) (1) 

where Sl r and <fi are the Rabi frequency and the initial phase of the pulse, respectively. Based on the Hamiltonian 
(1), it is straightforward to show that a pulse of duration t results in the following rotation 

|1) -> cosfi r i|l) -ie - ^sinfi r t|2), 

|2) -> -ie^ smQ. r t |1) + cosQ r t |2) . (2) 

Note that the resonant interaction can be done within a very short time by increasing the pulse Rabi frequency Q r 
(i.e., via increasing the pulse intensity). 

In the following, we also need the resonant interaction of the pulse with the |0) «-> |1) transiton of the qubit system 
1 [Fig. 1(c)]. In this case, we have 

|0) -> cosO r t|0) -ie"^sinf2 r i|l) , 

|1) -> -ie^sinft r t|0) + cosn r t|l), (3) 
where Q r is the Rabi frequency of the pulse. 

B. System-cavity resonant interaction 



Consider qubit system 1 coupled to a single-mode cavity field. Suppose that the cavity mode is resonant with 
the |1) <-> |2) transition while decoupled (highly detuned) from the |0) o |2) transition and the |0) <-> |1) transition 
of the qubit system l[Fig. 1(b)]. The interaction Hamiltonian in the interaction picture, after the rotating- wave 
approximation, is described by 

Ej = h(ga + |1) (2| +h.c.) , (4) 

where a + and a are the creation and annihilation operators of the cavity mode, and g is the coupling constant between 
the cavity mode and the |1) <H> |2) transition of the qubit system. Under the Hamiltonian (4), the time evolution of 
the states |2) |0) c and |1) |l) c of the whole system are as follows 

|2)|0) c -+ cos^|2)|0) c -isin^|l)|l) c , 

|l)|l) c -+ -i siiifft |2> |0> c + cosfft |1> |1> C , (5) 

where the states |0) c and |l) c are the cavity-mode vacuum state and single-photon state, respectively. 

C. System-cavity-pulse off-resonant Raman coupling 

Consider (n — 1) four-level qubit systems (2, 3, n). The cavity mode is coupled to the |1) -H- |3) transition of each 
qubit system, but decoupled (highly detuned) from the transition between any other two levels [Fig. 1(d)]. In addition, 
a classical pulse is applied to each one of qubit systems (2,3, ...,n), which is coupled to the |2) <H> |3) transition but 
decoupled from the transition between any other two levels [Fig. 1(d)]. In the interaction picture, the Hamiltonian 
for the whole system is 

n 

H = Y, ftaAe-^a+a^ + H.c.) + /<<>,;- ' ^s:,., + H.c.)] , (6) 

.7=2 

where the subscript j represents the jth qubit system, Vlj is the Rabi frequency of the pulse applied to the jth qubit 
system, cr^j • = |1),- (3| , ^23 j = (3| > and gj is the coupling constant between the cavity mode and the |1) o |3) 
transition of the jth qubit system. 

The detuning between the |2) -H- |3) transition frequency lj 32 of the jth qubit system and the frequency Uj of the 
pulse applied to the jth qubit system is Aj = lu 3 32 — ujj (j = 2, 3, n) [Fig. 1(d)]. In addition, the detuning between 

the |1) O 1 3) transition frequency oj j 31 of the jth qubit system and the cavity-mode frequency uj c is A c j = uj j 31 — u c 
[Fig. 1(d)]. Note that the detuning A c j is not the same for each qubit system (i.e., dependent of j) in the case when 
the level spacings are nonidentical for each qubit system. Under the condition A CJ 3> gj and Aj ftj, the level |3) 
can be adiabatically eliminated [26] and the effective Hamiltonian is thus given by [27-29] 
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where trf,, - |1>, (2| , X j - ^(l/A, + l/A cJ ), and 



&i — A CJ - - Aj — lo j 31 



(7) 



(8) 



From Eq. (8), it can be seen that the detuning Sj is adjustable by changing the pulse frequency ojj. With suitable 
choice of the pulse frequencies, we can satisfy 



5 n = 6, 



(9) 



which will apply below. 

For S 3> g 2 /A c .j, flj/Aj, \j, there is no energy exchange between the qubit systems and the cavity mode. Thus, 
under the condition (9), the effective Hamiltonian (7) can be written as [25,30,31] 

n /ft 2 - 2 



i=2 



^|2) J <2| + ^-o+a|l) i <l| 
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(a+a|l).(l|-aa+ 12)^.(21) 



(10) 



where the two terms in the second line above describe the photon-number dependent Stark shifts induced by the 
off-resonant Raman coupling, and the two terms in the last parentheses describe the "dipole" coupling between the 
two qubit systems (j,f) mediated by the cavity mode and the classical pulses. In the case when the level |2) of each 
qubit system is not populated, the Hamiltonian (10) reduces to 



j=2 



6 



a+a 



(11) 



It is easy to see that the states |0) ■ |0) c , |1) • |0) c , and |0) • |l) c remain unchanged under the Hamiltonian (11). However, 
if the cavity mode is initially in the photon state |l) c , the time evolution of the state |1) . of the jth qubit system 
under the Hamiltonian (11) is given by 



|1),|1> C 



A\j t 



1>,-|1>c 



(12) 



where A., = + ^f, which can be further written as 



1 A, 



c,3 



+ 2M, 

46 1 1 



(13) 



Note that the parameters gj is not adjustable once the qubit systems (e.g., solid-state devices) are designed and built 
in a cavity, and the detuning A c j = oj j 31 — lo c is fixed when the cavity mode frequency is chosen. As can be seen in 
Eq. (13), the parameter Xj here is adjustable by changing the pulse Rabi frequency Qj. 

III. GENERATION OF MULTI-QUBIT GHZ STATES 

Let us consider n qubit systems (1,2, ...,n) in a single- mode cavity. The qubit system 1 has three levels shown 
in Fig. l(a,b,c) while the four levels of the qubit systems (2,3, ...,n) are depicted in Fig. 1(d). For qubit system 1, 
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the cavity mode is resonant with the |1) -R- |2) transition but highly detuned from the transition between any two 
other levels [Fig. 1(b)]. In contrast, for qubit systems (2, 3, n), the cavity mode is off-resonant with the |1) O |3) 
transition but highly detuned from the transition between any other two levels [Fig. 1(d)]. These requirements can be 
achieved by an appropriate choice of qubit systems (e.g., atoms), prior adjustment of the level spacings of the qubit 
systems (e.g., superconducting devices), or prior adjustment of the cavity mode frequency before the operation. Note 
that the cavity mode frequency for both optical cavities and microwave cavities can be changed in various experiments 
(e.g., see, [32-36]). And, for superconducting qubit systems, the level spacings can be readily adjusted by varying the 
external parameters (e.g., the external magnetic flux and gate voltage for superconducting charge-qubit systems, the 
current bias or flux bias in the case of superconducting phasc-qubit systems and flux-qubit systems, see e.g. [37-39]). 

Suppose that the cavity mode and each of qubit systems (1, 2, n) are initially in |0) c and (|0) + |1)) /V2, respec- 
tively. The initial state for each of the qubit systems here can be easily prepared by the application of classical pulses. 
The whole procedure for preparing qubit systems (1, 2, n) in a GHZ state is shown as follows: 

Step (i): Apply a classical pulse (with a frequency uj — uj 2 i and (f> = — ir/2) to qubit system 1 for a duration 
tia = i"/ (2fi r ) [Fig. 1(a)], wait to have the cavity mode resonantly interacting with the |1) <-» |2) transition of the 
qubit system 1 for a time interval tu, = ■n/{2g) [Fig. 1(b)], and then apply a pulse (with a frequency uj = loiq and 
4> = — it/2) to qubit system 1 for a duration t\ c = n/(2il r ) [Fig. 1(c)]. According to Eqs. (2), (3), and (5), it can be 
seen that after the operation of this step, the following transformation is obtained: 



|0>! |0> c 
|l>i|0> £ 



|0) x |0> c 
|2>! |0> c 



|0>!|0> c 

-i|l>i|l>c 



|l)i|0) c 
i\0)i\l) c ' 



which leads the initial state Yl™ =1 + |0) c °f the whole system to the following state: 

n 

n (lO)^ + 11)^.) |0> c -h i jO)! |1> C ) . 



(14) 



(15) 



j=2 



Here and below, the normalization factor 1/2"/ 2 is omitted for simplicity. Eq. (15) shows that the levels |1) and |2) 
of qubit system 1 are not populated in the case when the cavity mode is in the single-photon state |l) c . Hence, after 
the operation of this step, the qubit system 1 is decoupled from the cavity mode during the operation of next step. 

Step (ii): Apply a classical pulse (with a duration t 2 ) to each of qubit systems (2, 3, n) to induce the off-resonant 
Raman coupling described in Sec. II C [Fig. 1(d)]. By adjusting the pulse frequencies, we set 5 2 = 83 = ... = 5 n = 5. 
Since the level |2) for each of qubit systems (2,3, ...,n) is not populated during the operation of step (i) above, the 
effective Hamiltonian describing this step of operation is given by Eq. (11). Accordingly, when the cavity mode is in 
the single-photon state |l) c , the time evolution of the state for qubit system j is then given by Eq. (12). Set 
\ 2 = A3 = • • • = A„ = A, which can be readily achieved by adjusting the Rabi frequencies fl 2 , ^3,---,^n of the 
pulses applied to qubit systems (2, 3, n). From Eq. (12), one can see that for t 2 — tt/X, we have the transformation 
|1) . |l) c - |1) . |l) c (j = 2,3, ...n). Thus, the state (15) becomes 



Tl Tl 

n (i°>, + 1 1 ),) i!>i i°)c + * n (i°>, - i 1 ), 



(16) 



j=2 



j=2 



which shows that the n qubit systems (1, 2, ...n) have been entangled to each other after the above operations. Since 
the qubit systems (1, 2, n) are also entangled with the cavity mode, we will need to perform the following operation 
to disentangle qubit systems (1, 2, ...n) from the cavity mode. 

Step (hi): Perform a reverse operation described in step (i). Namely, apply a classical pulse (with a frequency 

uj = u>io and an initial phase 4> — 7r/2) to qubit system 1 for a duration t\ c = ir/ (2fl T ^j [Fig. 1(c)], wait to have the 

cavity mode resonantly interacting with the |1) <H> |2) transition of the qubit system 1 for a time interval t\b — n/(2g) 
[Fig. 1(b)], and then apply a pulse (with a frequency oj — uj 21 and (f> — ir/2) to qubit system 1 for a duration 
tia = ft/(2Q r ) [Fig. 1(a)]. According to Eqs. (2), (3), and (5), it can be seen that after the operation of this step, we 
have: 



|l)i|0) c 

|0>! |1> C 



|0> 1 |0) c 

-|l>l|l> £ 



|0> 1 |0) c 
i|2>!|0) c 



|0) 1 |0) c 
i|l>i|0> £ 



which leads the state (16) to the following state 



n(K+iO ^i-nO^-i 1 ^)! 1 )! 

3=2 3 =2 



|o) c . 



(17) 



(18) 
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Eq. (18) demonstrates that the cavity mode returns to its original vacuum state and the qubit systems (1,2, ...,n) 
have been disentangled from the cavity mode after the above operations. 
The left part of the product in the state (18) can be rewritten as 

|0> 1 |+>2"-|+> n -|l>l|->2- ■■!->„. (19) 

where = [0)^. + |1) ., and [-^ = [0}^. - 1 1)^- (j = 2,3,...,n). Since is orthogonal to , the state (19) is a 
GHZ entangled state of n qubits. 

To reduce the operation errors, the level-spacing inhomogeneity in each four-level system needs to be larger than 
the bandwidth of the applied pulse, such that the overlapping of pulse spectra or the transition between any two 
irrelevant levels is negligible. This requirement can be achieved by prior adjustment of the qubit level spacings. For 
superconducting qubit systems, the level spacings can be readily adjusted by varying the external parameters [37-39]. 
To simplify our presentation, we will not give a detailed discussion here. Note that how to have the irrelevant levels 
not affected by the pulses or the cavity mode via prior adjustment of the level spacings was previously discussed (e.g., 
see [27,41]). 

Several additional points need to be made, which are as follows: 

(i) Since the cavity mode is off-resonant with the |1) <-> |3) transition of qubit system j (j = 2, 3, ...n) during steps 
(i) and (iii), a phase shift exp(iifj) happens to the state |1) of qubit system j when the cavity mode is in the single 
photon state |l) c for either of steps (i) and (iii). Here, ipj = g? (tu + t\ c ) /A CJ -. It is easy to find that if this unwanted 
phase shift for qubit system j is considered, the fidelity of the prepared GHZ state is given by 
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J=2 



'I + e-^i 



j=2 



(20) 



which shows that for ipj ~ 0, i.e., when the condition 

g) (tu + tu) /A cJ < 1 (21) 

is met, we have F <~ 1. The condition (21) can be reached by increasing g and f2 r to shorten the operation time 
tib + tie or increasing the ratio Acj/g?. 

(ii) For steps (i) and (iii), the pulse process depicted in Fig. 1(a) and (c) must be much faster than the process of 
the cavity mode resonantly interacting with the |1) <-> |2) transition of the qubit system 1 [Fig. 1(b)], i.e., 

tla,tlc«tl6, (22) 

such that the internal transition between the two levels |1) and |2) of qubit system 1 during the pulses, induced by 
the resonant interaction of the cavity mode with the |1) <-> |2) transition of qubit system 1, is negligible. Note that 
the condition (22) can be readily achieved by increasing the pulse Rabi frequencies fl r and fl r such that fi r , fl r ^> g. 

(iii) The level |2) of qubit system 1 is populated for a time interval t\ a + t u during step (i) or step (iii), thus the 
condition 

£la + £l6<72"r 1 >72 p 1 (23) 

needs to be satisfied in order to reduce decoherence caused due to spontaneous emission and dephasing of the level 
|2) of qubit system 1. Here, 7^. and 7^, are the energy relaxation time and dephasing time of the level |2) of qubit 
system 1, respectively. 

From the above description, one can see that this proposal has the following advantages: 

(i) No identical qubit-cavity coupling constants for the (n — 1) qubit systems (2,3, ...,n) are required, thus, this 
proposal is tolerant to the qubit-system parameter nonuniformity and nonexact placement of qubits in a cavity; 

(ii) No adjustment of the level spacings of the qubit systems or adjustment of the cavity mode frequency during 
the entire operation is needed; 

(iii) Neither auxiliary qubit systems nor measurement on the qubit states is needed; 

(iv) No photon detection is needed; 

(v) The entanglement preparation is deterministic; 

(vi) The entire operation time is given by 

t = n/g + n/Cl r + TT/tt r + 7r/A, (24) 



which is independent of the number of qubits in the cavity. 
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FIG. 2: (Color online) (a) System-cavity-pulse off-resonant Raman coupling for n qubit systems (1,2, ...,n) in a cavity. For 
simplicity, we here only draw a figure for qubit system j interacting with the cavity mode and a classical pulse (j — 1, 2, n). 
Sj — A c ,j — Aj is the detuning of the cavity mode with the pulse, A c ,j = uj 20 — u c is the detuning between the cavity mode 
frequency oj c and the |0) <-¥ |2) transition frequency oj 20 of qubit system j, Aj = uj 21 — ujj is the detuning between the pulse 
frequency ujj and the |1) |2) transition frequency ui 21 of qubit system j. The blue arrow dashed line represents a second pulse 
applied to the qubit system j to cancel the Stark shift of the level 1 1) induced by the first pulse, (b) System-cavity off-resonant 
interaction for each of n three-level qubit systems (1,2, ...,n) in Ref. [17]. g' is the non-resonant coupling strength between 
the cavity mode and the |0) <->■ |2) transition. 5 = UJ20 — w c is the detuning of the cavity mode frequency with the |0) o |2) 
transition frequency. 

In addition, for the description given above, it can be seen that: 

(vi) During the entire operation, the levels |2) and |3) of qubit systems (2,3, ...,n) are unpopulated and thus 
decoherence due to spontaneous emission and dephasing from these levels are greatly reduced; 

(vh) The level |2) of qubit system 1 is only populated for a very short time t\ a + tn, during step (i) or step (hi); 
and 

(viii) The level |2) of qubit system 1 is not occupied during step (ii) [note that the operation of this step (ii) requires 
much longer time than both step (i) and step (hi)]. 

Above we have discussed how to prepare a multi-qubit GHZ state with a three-level qubit system and (n — 1) 
four-level qubit systems in a cavity. The discussion given above is based on qubit systems for which the level spacings 
become narrower as the levels go up [see Fig. l(a,b,c) for qubit system 1 and Fig. 1(d) for qubit systems (2, 3, j)]. 
Note that this limitation is unnecessary. Namely, this proposal is also applicable to: (i) the qubit system 1 with such 
three levels, for which the level spacing between the two lowest levels |0) and |1) is smaller than that between the 
two upper levels |1) and |2) ; and (ii) the (n — 1) qubit systems (2, 3, n) with these four levels, for which the level 
spacing between the two levels \i) and \i + 1) is larger or smaller than that between the two levels \i + 1) and \i + 2) 
(here, i = 0, 1). Furthermore, since the level |0) of each of the qubit systems (2, 3, n) was not involved during the 
entire operation, one can choose the qubit systems (2, 3, n) for which the transition between the two lowest levels 
|0) and |1) is forbidden or weak to avoid or reduce decoherence caused by the spontaneous emission from the level 

|1>- 

The three-level or four-level qubit systems here are widely available in natural atoms and also in artificial atoms 
such as superconducting charge-qubit systems [37], phase-qubit systems [38,39], and fhix-qubit systems [37,40]. For a 
detailed discussion, see Ref. [27]. 

IV. COMPARING WITH PREVIOUS PROPOSALS 

In this section, we will give a comparison between our proposal and previous ones. For simplicity, we will compare 
our proposal with the ones in [17,25]. To the best of our knowledge, the proposals in [17,25] are most closely related 
to our work. 

A. Comparison with the proposal in [25] 

An n-qubit GHZ state can be prepared without real excitation of the cavity mode, by using the method introduced 
in [25]. To see this, consider n three-level qubit systems (1, 2, n) in a cavity. Assume that the cavity mode is coupled 
to the |0) O 1 2) transition and the pulses are coupled to the |1) O |2) transition, to establish the system-cavity-pulse 
off-resonant Raman coupling [Fig. 2(a)]. In this case, we can obtain an effective Hamiltonian, i.e., the Hamiltonian 
(10) with a replacement of |1) , |2) , and X)j=2 by |0)> 1-0 > anc ^ respectively. When the cavity mode is initially 
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in a vacuum state |0) c , this effective Hamiltonian reduces to 




(25) 

where Aj = uj^i — Wj, A c j = uj^o — w o and S = Sj = A c j — Aj (which can be reached via adjustment of the pulse 
frequencies). The notation of Xj 1S the same as that given in Sec. III. To prepare the qubit systems (1,2, ...,n) in a 
GHZ state, one will need to apply a second pulse to each of the qubit systems (1, 2, n) to cancel the Stark shifts, 
i.e., the first term of Eq. (25). After that, the Hamiltonian (25) becomes 

n 2 n 

(26) 

For the case of xi — X2 = ■■■ = Xn = X (achievable by changing the pulse Rabi frequencies), this Hamiltonian is the 
same as the Hamiltonian (6) presented in Ref. [25] for n two-level atoms interacting dispersively with a single-mode 
cavity. According to the discussion in Ref. [25], an n-qubit GHZ state can be prepared based on the Hamiltonian (26). 

From the description here, one can see that to achieve identical Raman transition strengths xi, X2, —, an d Xn, the 
pulses applied to different qubit systems should have different Rabi frequencies if the qubit-cavity coupling constants 
are nonidentical. And, to obtain identical detuning S for each qubit system, adjustment of the pulse frequencies would 
be needed. Furthermore, to generate GHZ states, a second off-resonant pulse should be applied to each qubit system 
to cancel the Stark shift induced by the first pulse. Namely, for generation of an n-qubit GHZ state, a total of 2n 
pulses would be required. In contrast, as shown above, our present proposal, which employs four-level qubit systems, 
needs only (n — 1) pulses applied to the qubit systems (2, 3, n), i.e., one pulse for each of qubit systems (2, 3, n). 
Hence, the use of the pulses is significantly reduced in the present proposal. 

B. Comparison with the proposal in [17] 

In Ref. [17], an auxiliary qubit system a was used to create a single photon in the cavity mode, which was then 
employed to prepare the n qubit systems (1, 2, n) in a GHZ state. As discussed there, measurement on the states 
of the auxiliary qubit system a and adjustment of the level spacings of each qubit system were both needed during 
the entire operation. In addition, qubit systems (1,2, ...,n) were required to have identical three- level structures and 
the qubit-cavity coupling constants for qubit systems (1, 2, n) were needed to be the same. 

The cavity mode in [17] was set to be off-resonant resonant with the |0) o |2) transition of each of qubit systems 
(1, 2, n), with a detuning S = W20 — w c [Fig. 2(b)]. The coupling constant between the cavity mode and the 
|0) <-» 1 2) transition for each of qubit systems (1,2, ...,n) was denoted as g' . As shown in [17], a multiqubit GHZ 
state was generated by simultaneously performing a common phase shift e™ on the state |0) of each of qubit systems 
(1,2, ...,n) with assistance of the cavity photon, by having e tXt — e m (i.e., Xt = n). Here, A = (g') 2 /S, which was 
originally defined in [17]. For a detailed discussion, see Ref. [17]. 

It is noted that the method in [17] can not work in the case when the qubit-cavity coupling constants are nonidentical. 
The reason is that a common phase shift e l7T can not be simultaneously performed on the state |0) of each of qubit 
systems (1, 2, n). For instances, when the qubit-cavity coupling constants g'j and g' k for qubits j and k are not the 
same (resulting in different Aj — {g'j) 2 /5 and A^ = (g' k ) 2 /5), one can not have both of Xjt and Xkt to be equal to tt 
for a given time t. The discussion here applies to identical atoms in a cavity, for which the detuning 5 of each atom 
with the cavity mode is the same but the qubit-cavity couplings are nonidentical due to nonexact placement of atoms 
in the cavity. 

Let us now see if the method in [17] can work for a situation that the detuning of each qubit system with the 
cavity mode and the qubit-cavity coupling constants are both nonidentical. This situation applies to superconducting 
qubit systems in a cavity. As discussed above, to prepare a multiqubit GHZ state, the condition Xj = needs to 
be satisfied for arbitrary two superconducting qubit systems j and k. Here, Aj = (g'j) 2 /5j and A^ = (g' k ) 2 /^k- Note 
that once the superconducting qubit systems are designed and the cavity mode frequency is chosen, the detunings Sj 
and 5k are fixed. Thus, in order to obtain Aj = Afe, i.e., (g'j) 2 /Sj — (g' k ) 2 /Sk, it would be required to exactly place the 
qubit systems in the cavity to have desired qubit-cavity coupling constants g'j and g' k , which is not easy to achieve in 
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experiments due to the fabrication error. This problem becomes more apparent especially when the number of qubit 
systems in a cavity is large. In contrast, as shown above, this difficulty is avoided in our present proposal, because the 
effective coupling strength Aj in Eq. (13) can be adjusted by changing the intensity of the pulses, instead of exactly 
placing qubit systems in a cavity to have desired qubit-cavity coupling constants. 

From the discussion here, it can be concluded that our present proposal is quite different from the previous ones 
[17,25]. Due to the use of a four-level structure, the present proposal can be used to prepare a multiqubit GHZ state 
for nonidentical qubit-cavity coupling constants (compared with the proposal in [17]), and requires less application of 
pulses (compared with the proposal in [25]). 

V. DISCUSSION 

In this section, we will give a brief discussion on the experimental issues. For the method to work: 

(i) The occupation probability pj of the level |3) for qubit system j (j = 2, 3, n) during step (ii) is given by [28] 

2 2 
Pj ~ 4+(A J /fi J ) 2 + 4+(A c , j /< 7j ) 2 ' 

which needs to be negligibly small in order to reduce the operation error. 

(ii) As discussed above, the conditions (21), (22), and (23) need to be satisfied; 

(iii) The total operation time r given in Eq. (24) should be much shorter than the energy relaxation time 7 1 ^ 1 and 
dephasing time 7-jl 1 of the level |1) and the lifetime of the cavity mode = Q/2ttv c , where Q is the (loaded) quality 
factor of the cavity. 

The above requirements can in principle be realized, since one can: (i) reduce pj by increasing the ratio of Aj/Qj 
and A c j/gj; (ii) shorten ti a ,ii6,iic by increasing the pulse Rabi frequency fi r , Cl r , and the coupling constant g; (iii) 
design qubit systems (e.g., superconducting devices) to have sufficiently long energy relaxation times 7^. and 7^. and 
dephasing times 7^ and 7^ or choose qubit systems (e.g., atoms) with long decoherence times of the levels |1) and 
1 2) ; (iii) increase k _1 by employing a high-Q cavity so that the cavity dissipation is negligible during the operation. 

For the sake of definitiveness, let us consider the experimental possibility of generation of a six-qubit GHZ state with 
six superconducting qubits coupled to a resonator (Fig. 3). With the choice of fl r , fl r ~ lOg, we have t\ a , t\ c <~ 0.057r/g 
and tib ~ 0.5n/g. By setting Aj ~ lOilj, A C J ~ Wgj, and fl, ~ 0.9<7j, we have 5 = 5j ~ 10g^ / A c j <~ 10f2j/Aj <~ 10;\;j, 
leading to Aj <~ O.llgj. For gj <~ 0.2g, we can set A 2 = A 3 = ... = = A ~ 0.022g, which can be achieved by adjusting 
the Rabi frequencies Q2, ^3, ^6 of the pulses applied to the qubit systems (2, 3, 6) as discussed above. Thus, we 
have t ~ 46. 671-/(7. As a rough estimate, assume g/2n <~ 220 MHz, which could be reached for a superconducting qubit 
system coupled to a one-dimensional standing- wave CPW (coplanar waveguide) transmission line resonator [42]. For 
the g chosen here, we have: (i) t lai t lc <~ 0.1 ns and tu ~ 1.1 ns, which shows that the condition (22) is satisfied; (ii) 
ha + hb ~ 1-2 ns, showing that the condition (23) is met for min{7 2r 1 , 7^} ~ 0.5 fxs [38,43]; and (iii) r ~ 0.1 (is, 
much shorter than min{7^ r 1 , 7^} ~ 1 (is [38,43]. In addition, consider a resonator with frequency v c ~ 3 GHz (e.g., 
Rcf. [44]) and Q ~ 5 x 10 4 , we have kT 1 ~ 2.5 /is, which is much longer than the total operation time r. Note that 
superconducting coplanar waveguide resonators with a quality factor Q > 10 6 have been experimentally demonstrated 
[45]. 

For the choice of Aj ~ lOfij and A CJ - ~ 10gj here, we have pj ~ 0.04, which can be further reduced by increasing the 
ratio of Aj/flj, and A c j/gj. In addition, for the parameters chosen here, when the unwanted phase shifts mentioned 
above are considered, a simple calculation shows that the fidelity F is <~ 0.992. We should mention that further 
investigation is needed for each particular experimental set-up. However, this requires a rather lengthy and complex 
analysis, which is beyond the scope of this theoretical work. 

As shown in Sec. Ill, this proposal requires simultaneous application of pulses with different frequencies and intensi- 
ties during the operation of step (ii). For the set up in Fig. 3(a), simultaneous application of several pulses of different 
frequencies and intensities to qubit systems 2,3,4, and 5 is feasible in experiments [46]. According to Ref. [46], for a 
qubit system containing a superconducting loop, one can place an AC current loop on the qubit superconducting loop 
to create an AC flux $ e (i.e., a classical magnetic pulse) threading the qubit loop [Fig. 3(b,c,d)]; and the frequencies 
and intensities of pulses applied to qubit systems 2, 3, 4, and 5 can be readily and simultaneously changed by varying 
the frequencies and intensities of the AC loop currents of the qubit systems 2, 3, 4, and 5 at the same time. 

VI. CONCLUSION 

In summary, we have proposed a way for creating n-qubit GHZ entangled states with a three-level qubit system 
and (n — 1) four-level qubit systems in a cavity or coupled to a resonator. The main advantage of this proposal is 
that no identical qubit-cavity coupling constants are required. As a result, qubit systems (e.g., solid-state qubits), 
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Ej aEj 
(b) (c) (d) 

FIG. 3: (Color online) (a) Sketch of the setup for six superconducting qubit systems and a (grey) standing-wave quasi-one- 
dimensional coplanar waveguide resonator. The two blue curved lines represent the standing wave magnetic field, which is in 
the z-direction. The green dot represents the qubit system 1, which is placed at an antinode of the standing wave magnetic 
field to achieve the maximal qubit-cavity coupling constant g. The red dots represent the qubit systems 2, 3, and n, which 
are placed at locations where the magnetic fields are the same. Each qubit system could be a superconducting charge-qubit 
system shown in (b), flux-biased phase-qubit system in (c), and flux-qubit system in (d). The superconducting loop of each 
qubit system, which is a large square for (b) and (d) while a large circle for (c), is located in the plane of the resonator between 
the two lateral ground planes (i.e., the x-y plane). Each qubit system is coupled to the cavity mode through a magnetic flux 
$ c threading the superconducting loop, which is created by the magnetic field of the cavity mode. A classical magnetic pulse 
is applied to each qubit system through an AC flux <3? e threading the qubit superconducting loop, which is created by an AC 
current loop (i.e., the red dashed-line loop) placed on the qubit loop. The frequency and intensity of the pulse can be adjusted 
by changing the frequency and intensity of the AC loop current. Ej is the Josephson junction energy (0.6 < a < 0.8) and V g 
is the gate voltage. A is the wavelength of the resonator mode, and L is the length of the resonator. 

which often have parameter nonuniformity, can be used; and no exact placement of qubits in a cavity is needed. 
Note that for solid-state qubits (e.g., superconducting qubits), it is difficult to design qubits with identical device 
parameters; and it is experimentally challenging to place many qubits at different locations where the magnetic fields 
or electric fields are exactly the same. These hardships are avoided in this proposal. Another interesting property of 
this proposal is that during the entire operation, there is no need of adjusting the level spacings of the qubit systems 
or the cavity mode frequency. In experiments, adjustment of the qubit level spacings or the cavity mode frequency 
during the operation is not desired; and extra errors may be induced by adjustment of the qubit level spacings or the 
cavity mode frequency. Furthermore, as shown above, this proposal has these additional advantages: (i) The n-qubit 
GHZ state is prepared deterministically and the operation time is independent of the number of qubits in the cavity; 
(ii) Neither auxiliary qubit systems nor measurement on the qubit states is needed; and (iii) No photon detection is 
needed. This proposal is quite general, which can be applied to various types of superconducting qubits and atoms 
trapped in a cavity. 
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